Rashba coupling in quantum dots: exact solution 
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We present an analytic solution to the problem of the Rashba spin-orbit coupling in semiconductor 
quantum dots. We calculate the exact energy spectrum, wave-functions, and spin-flip relaxation 
times. We discuss various effects inaccessible via perturbation theory. In particular, we find that 
the effective gyromagnetic ratio is strongly suppressed by the spin-orbit coupling. The spin-flip 
relaxation rate has a maximum as a function of the spin-orbit coupling and is therefore suppressed 
in both the weak- and strong coupling limits. 

PACS numbers: 71.70.Ej, 73. 21. La, 72.25.Rb 



^ I In recent years, spin-orbit (SO) effects in quantum dots attracted much attention, as it has become clear that these 
effects play a crucial role in the novel field of spintronics?^ . 

The combination of a confined geometry and SO coupling has interesting consequences for the electron spectrum^. 
. Since future quantum computation devices would have to control coherent spin states over sufficiently long time- 
scales"^, it is important to understand spin relaxation mechanisms, most of which are rooted in the SO coupling. 

SO interactions can arise in quantum dots by various mechanisms related to electron confinement and symmetry 
breaking and are generally known as the Rashba term^- and the Dresselhaus term^. For most experimental realisations, 
'"^ [ quantum dots can be described as effectively two-dimensional systems in a confining potential which is usually 
^ ' modelled as hard- wall or harmonic confinement. 

O ] To our knowledge all existing theoretical studies of SO effects in such systems rely on perturbative schemes or 
numerical simulations. The purpose of this paper is to provide an exact solution of the quantum mechanical problem 
of the combined effects of the SO coupling, confinement, and magnetic field. 
Cn| . In the bulk (i.e. without confinement), the problem was solved in the original paper—, see also&, while the effects of 
K*" ' the confining potential (but without SO coupling) were studied in Ref—. We shall show here how these exact solutions 
\l , can be combined and generalised. 

' The one-particle Hamiltonian describing an electron in a two-dimensional quantum dot is of the form: 

O ■ p2 1 

^1 H = — + V{p) + ajiip^ay - Py(Tx) + -gfisBaz , (1) 

■ where m is the effective electron mass, g is the effective gyromagnetic ratio, and /i^ is the Bohr magneton. A constant 
magnetic field B is introduced via the Zeeman term above and the Peierls substitution, p = — iV — |A (we use the 
axial gauge, Ap = and A^p = Bp/2); an is the strength of the SO coupling, the Pauli matrices are defined as 
standard, the electron charge is e = — |e|, and we set h = 1. In V{p) is the (symmetric) confining potential. In 
this paper we will consider a hard-wall confining potential, i.e. V{p) = for p < i? and V{p) = oo for p > R, R being 
^ ■ the radius of the dot. We have included the Rashba term rather than the Dresselhaus term, which would be of the 
Q form aoiPxO'x ~Pv<^y)- (The Rashba term maybe dominant, since the coupling strength an can be varied by system 
O ■ design, e.g., values of aji ~ 6 x 10~^eV cm were reported for InAlAs/InGaAs structures-, whereas the typical values 
^1 of au are about IO~^eV cm^.) These terms transform into each other under the spin rotation: cf^ ^ (Jy, Oz —cfz- 
So our results will only need a trivial modification in the case when a solo Dresselhaus term is present. 

Hamiltonian commutes with the z-projection of the total momentum operator = h + \'^z\ Iz = ^id^ 
(assuming the axial gauge). The operator jz is therefore conserved. The eigcnfunctions of the total momentum 
operator, with a half-integer eigenvalue j, are of the form: 

^,(p, ^) = {e^^^-'/'^'^Mp), e'(^"+i/2)'^.9,(p)). (2) 

In zero field there is an additional symmetry jz — > —jz related to time inversion. The states with the projections of 
momenta equal to j and —j are Kramers doublets. Introducing the operators, Vj^j = ±^ ~ ji ~ ' 



^ pctpYdp) p^y 2c ^ 



the Schrodinger equation becomes. 
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^T-i/2fi + 2m - y^JiBB - V{p)\ /, - 2ma^V^^^.^,/2 5, = , 
^f+\/293 + 2™ + \91^bB - V{p)] g, - 2maflV^^^._,/2 /, = . 



(3) 



We first assume that the magnetic length as = \/c/\e\B is large compared to R, so that the orbital effects of the 
magnetic field can be neglected. Working with the dimcnsionless coordinate x = p/R, the system of equations (0) 
becomes 

+ e - h)f, - /3«vL°i+i/2<?, = , 
(A;+i/2 + ^ + - /3«v(^"i_,/,/, = , 

supplemented by the boundary conditions fj{l) = = 0. We have introduced two dimcnsionless parameters, 

= 2ajimR and h = mgpgR^B, characterising the strength of the SO coupling and the Zeeman term, respectively. 
The energy parameter is e = 2mER^. 

In the absence of the Rashba term and confinement potential, the solutions regular at the origin are fj{x) ~ 
Jj-i/2{kx) and gj{x) ~ Jj^i/2{kx) with fc^ = e ± /i, where Ji{x) are the Bessel functions. Due to the standard 
recurrence relations 

+ Jj±i/2{kx) = kJj^i/2{kx) , (5) 

the Rashba term simply acts as rising or lowering operator on the Bessel functions basis^. Therefore the following 
ansatz 



{fj{x)^9j{x)) = {diJj-i/2{kx),d2Jj+i/2{kx)) (6) 

solves the bulk problem in the presence of the SO coupling, provided that the coefficients di^2 satisfy the eigenvalue 
equation: 

' k'^ - e + h -l3Rk 
—Pnk k^ — e — h 

When considering the electron confined to the disk, it is seemingly impossible to impose the vanishing boundary 
conditions on the ansatz ® as Bessel functions with different indices are involved. Note, however, that as long as 
either I3r or h is non-zero, the bulk spectrum has two branches: e = k'^ ± \/0jJt^\-h^ . Therefore for a given value 
of e there are, in fact, two non-trivial solutions for the momentum k, 

,2 (2e + /3|)±v//3j; + 4e/3^ + 4/t2 
k±- ^ . 

We choose the amplitude ratios as /d^ = a+ and d2 /d^ = a_, where a±(e, (}R,h) — (3Rk±/{kj^ — €±h) . We are 
now able to satisfy the boundary conditions by combining these two linearly independent degenerate solutions: 

(/i(a;),5j(a;)) = d+{a+Jj_i/2{k+x),Jj+i/2{k+x)) 

+d~{Jj_i/2{k.x),a-Jj+i/2{k-x)) . (8) 
The boundary conditions lead to the eigenvalue equation 

FJj+i/2{k+)Jj^i/2{k-) + Jj+i/2{k-)Jj-i/2{k+) = , 

where the function F = — 0^k+k^ / {{k'^ — e + h)(k'^ — e — h)) . Notice that in the energy region |e| < h, fc_ becomes 
purely imaginary k- — > ik-, and the respective Bessel function becomes a modified one, I{k-x). 

In the absence of the Zeeman term we have k± = + /?^/4 ± /3_r/2 and df = ±df , so that F = 1 in the above 
equation. In this case, the equation is invariant under the change j — > —j (Kramers degeneracy). At /3_r = all states 
with I ^ are four-fold degenerate, while I = states are doubly degenerate. According to the standard analysis, 
the SO coupling splits all the I ^ states into two Kramers doublets with j = Z -I- 1/2 and j — I — 1/2, while I = 
states naturally remain Kramers doublets. The specifics of the Rashba term is that, at small (3r, the SO splitting 
starts at the order j3jf . The evolution of the first few energy levels with the parameter Pr is shown in Fig.^ (Values 



di 
d2 



= 0. 



(7) 



3 



0.2 - 




LU a — ^ 

I , \ , \ , I , L 

0.5 1 1.5 2 



Spin-orbit coupling 

FIG. 1: Rashba splittings. Energy as a function of Pr for the states (1/2, 0) (a), (3/2, 0) (b), (1/2, 1) (c), (5/2, 0) (d), (3/2, 1) 
(e), and (1/2, 2) (f). {R = lOnm, B = 0.) 
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FIG. 2: Rashba splittings. Energy differences E{l/2, 1) - S(3/2, 0) and -E(3/2, 1) - -E(5/2, 0) as a function of I3r. Dashed lines 
correspond to the same quantities calculated at the second-order of perturbation theory. [R = lOnm, B = 0.) 



of m « 0.05me and g k, —0.4, typical for A3B5 structures, are used.) We label the energy eigenstates as (j, n) where 
n is a non-negative integer such that Ej^n < Ej^n+i at — 0. 

A comparison of the exact splittings to the perturbative ones is shown in Fig. |21 As one can see from this figure, 
the perturbation theory seriously (by 20-30%) overestimates Rashba splittings for realistic values of the parameter 

Upon inclusion of the Zeeman term, all Kramers doublets are also split so that all the degeneracy is completely 
lifted. Because of inherently small values of the gyromagnetic ratio g in the semiconductor quantum dots, for realistic 
magnetic fields the Zeeman splittings are small (10~^-10~^ meV ) in comparison to the characteristic energy separation 
between the levels. Because of the same reason, all the Zeeman splittings can be regarded as linear in the magnetic 
field. So, for the Zeeman splitting of the j's eigenstate we may write Sej — 2hfj{Pii) , where the function fjiPn) 
[/(O) = 1] pl^^ys the role of an effective gyromagnetic ratio, which non-trivially depends on Pfj. Indeed, expanding 
the eigenvalue equation in h, we find the following analytic formula for the gyromagnetic factor: 

where the prefactor is i = 4(e -I- /Sj^/A)/ {/Sue); e and k± are zero-field solutions. Numerical calculations show that 
/(/3-r) is suppressed by about 50% when Pn reaches 1.5. Increasing Pji suppresses the Zeeman splitting because the 
SO coupling entangles the spin degree of freedom, making it more difficult to polarise. 

The SO coupling is the main intrinsic mechanism for electron spin-flip transitions in quantum dots^i. In previous 
calculations of the spin-fiip rates, the SO coupling was considered as a perturbation, so that the electron spin and 
angular momentum were assumed to be independently conserved. In the full theory this is not the case. The 
'spin-flip' transitions in fact occur between the states j and —j. No such transition is possible within a degenerate 
Kramers doublet (Van Vleck cancellation). In the external magnetic field, the states j and —j are split by the 
Zeeman interaction. The SO coupling allows then for phonon assisted transitions between the Zeeman sublevels (of 
a given Kramers doublet). We concentrate on the most interesting case and consider the spin-flip transition rate 
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FIG. 3: The ground-state spin- flip relaxation rate W as a function of /3r for B = 1.5T and R — 10 ran. 



between the Zeeman sublevels of the ground state (j = ±1/2). Acoustic phonons dominate these processes at low 
temperatures. We consider piezoelectric interaction between the electrons and the acoustic phonons. The rate for the 
one-phonon transition within Zeeman sublevels is given by the Fermi golden rule. In analogy Ref.^, we obtain (at 
zero temperature): 

fe hu rY (AEzf 
W = ^ K, 



F= / cos4> ip*_^,2{x,(l)) il^ii2{x,(l))x^dx d4> (9) 

Ja Jo 

where = 3/2sf -|- 1/sf , K — 8/1057r, AEz is the energy difference between the states involved, hi4 is the non-zero 
component of the piezo-tensor and po is the mass density. The exact evolution of the spin-flip transition rate with 
the parameter for fixed h is shown in Fig. |31 We use typical material parameters for GaAs-type structures^". 
An interesting feature of the exact solution is the emergence of a maximum in the transition rate as a function of 
the spin-orbit coupling. The physical explanation is that while the Zeeman energy splitting decreases with (3ji, the 
electron-phonon matrix element saturates. For small Pr and h, we obtain F P^h so that the transition rate is 
W ~ h^0R, in accordance with the perturbative result of RefiS. 

So far we have assumed that as > R, which is the case for most experimental set-ups. For bigger dots and stronger 
magnetic fields, however, the orbital effects of the magnetic field can not be neglected. Fortunately, because of the 
very nature of the Peierls substitution, which has to be performed both in the kinetic energy term and in the SO 
term, the above analytic solution can be generalised to this case. 

Using the dimensionless variable ^ = p^/{2a\) we propose the ansatz 

- (di$(en,j,0,rf2<i>(eo + l,J + l,0) , 

where 

'i>(eo,J,0 = c,C^Af(-eo + ^^,|/| + l,C) 

^ / r(eo-Hn + i) \ ^-i)-^ 
l^r(eo-^ + i); r(|;| + i) 

with I = i — \ and M(a, c;^) is the confluent hypergeometric function (see^^). The factors Cj are inspired by the 
normalisation in the standard Landau problem and Cq is the 'energy' parameter to be determined. As before, we find 
two non-trivial solutions for it: 

± _ 2e - 1 + 7^ ± ^/(t^ + lY + 4572 + 4(s2 - s) 
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Magnetic field (T) 

FIG. 4: The spectrum in the orbital field. Energy as a function of magnetic field for the states (1/2,0) (al), (3/2,0) (bl), 
(1/2, 1) (cl), (5/2, 0) (dl), (3/2, 1) (el), and (1/2, 2) (fl), (-1/2, 0) (a2), (-3/2, 0) (b2), (-1/2, 1) (c2), (-5/2, 0) (d2), (-3/2, 1) 
(e2) and (-1/2, 2) (f2). (R = 50nm, f3R = 3.35.) 



We have introduced dimcnsionless parameters: e — E/lUc — 1/2 [uJc = \e\B/mc) for the energy (not to be confused 
with the electron charge), 7 — ad{2m/wc)^/'^ for the SO coupling, and s — gra/{Ame) for the Zeeman coupUng. Here 
rUe is the electron mass. The boundary conditions then lead to the eigenvalue equation: 

^(eo,^,^o)^(eo + + ^0) ^ (e" + 1)^ e+ - e + s 
$(e^,Z,Co)$(eJ + !,/ + !, Co) (e^J" + 1)* -e + s' 

where = ^^/(2a|.). This equation provides all the information about the energy spectrum of the system. We 
have investigated it numerically. The results are shown in Fig. ^ Note that the Zeeman splittings are still small for 
realistic fields. The Kramers doublets therefore survive the orbital field as long as there is no SO coupling. It is the 
combined effect of the orbital field and the Rashba term that lifts the Kramers degeneracy. 

To conclude, we presented an analytic solution to the problem of an electron in a quantum dot in the presence 
of both the magnetic field and SO coupling. We calculated various quantities of physical interest. For realistic 
parameters, the Rashba energy splittings are overestimated in perturbation theory. There is also a strong suppression 
of the effective gyromagnetic ratio by the SO coupling. The spin-flip relaxation rate has a maximum as a function 
of the SO coupling, a prediction that would be interesting to verify experimentally. Inclusion of the orbital magnetic 
field gives rise to a rich magneto-optical spectrum. We hope that our solution can be used in future research for 
obtaining further interesting results on the SO effects in quantum dots. 

We are grateful to Levitov, who has independently arrived at a similar solution with Rashba^^, for interesting 
discussions. G.L.'s and A.O.G.'s research is supported by the EPSRC grants GR/N19359 and GR/R70309 and the 
EU training network DIENOW. E.T.'s research is supported by the Centre for Functional Nanostructures of the 
Deutsche Forschungsgemeinschaft within project A2. 

Note added After completion of this work we learnt that the particular case of zero magnetic field was analised by 
Bulgakov and Sadree\*i^. 
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